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. Abstract. Given hydric capacity and nutrient now of a chemostat-like system, we analyse the influence 

\Q ' of a spatial structure on the output concentrations at steady-state. Three configurations are compared: 

perfectly-mixed, serial and parallel with diffusion rate. We show the existence of a threshold on the input 
concentration of nutrient for which the benefits of the serial and parallel configurations over the perfectly- 
mixed one are reversed. In addition, we show that the dependency of the output concentrations on the 
diffusion rate can be non-monotonic, and give precise conditions for the diffusion effect to be advantageous. 
The study encompasses dead-zone models. 
' Key-words. Chemostat model, interconnection, diffusion, global stability. 
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1 Introduction 



The chemostat is a popular apparatus, invented simultaneously by Monod [5D] and Novick & Szilard |23j . 
for the so-called continuous culture of micro-organisms. It has the advantage to study bacteria growth at 
steady state, in contrast to batch cultivation. In the classical experiments, the medium is assumed to be 
(■ — , , perfectly mixed, that justifies mathematical models described by systems of ordinary differential equations 

£f) ■ [29]. The chemostat model is also used in ecology for studying populations of micro-organisms, such as 

lake plankton or wetlands ecosystems. In natural ecosytems, or in industrial applications that use large 
bioreactors, the assumption of perfectly mixed medium is questionable. This is why spatial considerations 
have been introduced in the classical model of the chemostat, such as the gradostat model [16] that is a series 
of interconnected chemostats (of identicalvolumes). Segregated habitats are also considered in lakes, where 
the bottom can be modeled as a dead zone and nutrient mixing between the two zones is achieved by diffusion 
rate [3T] . The consideration of dead zones is also often used in bioprocesses modelling [TSl HH [SJ 123 [32 H7J ■ 
Series of chemostats, instead of single chemostat, have shown to potentially improve the performances 
of bioprocesses, reducing the total residence time [13j [TT1 [101 EH HI] or allowing species persistence [30J [24] . 
These properties have of course economical impacts for the biotechnological industry, and there is a significant 
literature on the design of series of reactors and comparison with plug-flow reactors (that can be seen as 
the limiting case of an arbitrary large number of tanks of arbitrary small volumes) [32[ [TJ \22\ [2J El [4j [5] . 
Sometimes a radial diffusion is also considered in plug- flow reactors [7J, but surprisingly, configurations of 
tanks in parallel have been much less investigated |15) . One can argue that knowing input rates and volumes 
of tanks in parallel, there dynamical characteristics can be studied separately, and there is no need of devoting 
a specific study for these configurations. This is no longer the case if one considers a passive communication 
between the tanks, through a membrane for instance. In saturated soils or wetlands, a spatial structure could 
be simply represented by separated domains with diffusive communication. This consideration is similar to 
patches models or islands models, commonly used in ecology [TBI [9] , or lattice differential equations [28] . 
For instance, a recent investigation studies the influence of such structures on a consumer/resource model 
[8]. Consumer/resource models in ecology are similar to chemostat models, apart the source terms that are 
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modeled as constant intakes of nutrient, instead of dilution rates that one rather met in liquid media. 

In this paper, we propose to bring new insight on parallel configurations of chemostats with communi- 
cation, in a spirit different than the one usually taken in bioprocesses design. One usually chooses a target 
for the output concentration of substrate, and looks for minimizing the total volume, or equivalently the 
residence time, among all the configurations that provide the same desired output at steady state. Here, we 
fix both the total hydric volume and the input flow and study the input-output map at steady-state, investi- 
gating the role of the spatial structure on the performances of the system. The performance is here measured 
by the level of substrate that is degraded by the system, and collected at the output. We draw precise com- 
parisons between the three configurations: perfectly mixed, serial and parallel (with diffusion rate) with the 
same total hydric volume and flow rate. This set of configurations is far to be exhaustive, being limited to 
two compartments only, but it is a first attempt to grasp this input-output map of a structured chemostat, 
and study how a spatial structure can modify this map, and what are the key parameters. We believe that 
this study is of interest for the modelling of ecosystems such as saturated soils for which it is not easy to know 
the spatial structure, and where one has only access to input-output observations of the substrate degradation. 

The paper is organized as follows. In Section [21 we present the three configurations under investigation 
and give the equations of the models. The main part of the paper is devoted to the analysis of the steady 
states, given in Section [3l The proofs of the global stability of the equilibriums are postponed to the 
Appendix, for lightening the presentation. Finally, discussion and numerical simulations are given in Section 
® 

2 The models 

The flow rate is labeled Q and V is the total capacity of the system. The three simple patterns we analyze 
are depicted on Figure Q] 

- one single compartment of volume V 

- two compartments of volume Vi , V2 such that V = V\ + V2 connected in cascade 

- two compartments of volume V\, V2 such that V = V± + V2 in parallel with a diffusive connection. 




one compartment serial connection parralel connection 



Figure 1: The set of configurations under investigation. 
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We recall the dynamical equations of resource (nutrient) and biomass concentrations, respectively denoted 
by Si and Xj in a compartment i of volume V, fed from a compartment z_ with a flow rate Qi and connected 
by diffusion rate d to a compartment id (see Figure [5]) . 

Si = " — ~Xi + -rr(Si_ — Si) + —{Si, — Si) 

y V V 

Vi Vi 




Figure 2: Possible interconnections of a compartment. 



For sake of simplicity of the analytical analysis, we assume that the growth function /i(-) is a linear function 
of the resource concentration: 

n(S) = mS 

In Section [4j we shall consider Monod growth function and show that the qualitative results of our study 
are not changed. The yield coefficient y of the bio-conversion is kept equal to one (this is always possible by 
choosing the unit measuring the biomass). It is convenient to write dimensionless concentrations: for each 
concentration Ci in the compartment i (Ci can denote Si or Xi), we define 

V r 

Ci — 771 — L/j , 

and 

= Yi 
n v ' 

We shall also consider that the time t is measured in units such that Q — V . Finally, we assume that the 
input concentration Si n is large enough to avoid the (trivial) wash-out equilibrium to be the only steady-state 
in each compartment. 



3 Steady-state analysis of the three configurations 

3.1 Configuration with one compartment 

The dynamical equations of the configuration with a single compartment are 

I S — SX ~\~ Sifi S 

\ X = SX — X 

The non-trivial equilibrium is (I, Sj„ — 1) under the condition Si n > 1. Then, one has 

-k i 

s out — 1 • 

Remark. This is a well known property from the theory of the chemostat that the output concentration at 
steady state is independent of the input concentration, provided this latter to be large enough (i.e. Si n > I). 
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3.2 Serial connection of two compartments 

The dynamical equations of the model with two compartments in series (see Figure [T]), assuming r to be 
different to and 1, are 

-SiXi + -(S in - Si) 



Si 

±1 

i-2 



1 



SiXl - 
-S2X2 
S2X2 + 



-Xl 



1 



■(si - s 2 ) 



(1) 



1 -r 

(xi-x 2 ) 



1 - r 



with r = Vi/V. 



Proposition 1. When Sj„ > 1/r, there exists an unique equilibrium (s*, x*, sj, x£) of (fT]) on the positive 
orthant. One has necessarily s* = 1/r and < min(l/r, 1/(1 — r)). Furthermore, one has 



< 1 



> 1 + 1/r 



Proof. One can readily check that there exists a non-trivial equilibrium (1/r, Sj n — 1/r) for the first 
compartment exactly when Si n > 1/r. Furthermore, this equilibrium is unique. Then, any equilibrium for 
the overall system ([T]) has to be {s\, Sm — s£) f° r the second compartment, with s£ solution of the equation 



S 2 {Sin - s 2 ) 



1 -r 



(1/r - s 2 ) 



(2) 



with S2 < l/ r - One can easily verify that there exists a unique solution of © on (0, 1/r). Graphically, 
S2 is the abscissa of the intersection of the graphs (see Figure |3]) of the polynomial function 



and the affine function 



iM = 



82) = s 2 (s 

1 



1-r 



(1/r -s 2 ) 



l/r-s 2 
















\ s 2 (s in -s 2 ) 








s 2 


1/r 


Si " S 2 



Figure 3: Graphical determination of sj- 

Remark that s; rl > 1/r implies the inequality <^(1/(1 — r)) > Z ( 1/ ( 1 — r)), from which one deduces s 2 < 
1/(1 — r). Finally one can compare s ou t = s 2 with the value obtained in the configuration of one compartment: 



S*out < 1 



(1) > 



> 1 + 1/r 



□ 
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3.3 Parallel interconnection of two compartments 

The dynamical equations of the model with two compartments in parallel and diffusion (see Figure [1} , 
assuming r to be different to and 1, are the following 

h = -six 1 + -(s in -sx) + -(s 2 -s 1 ) 
r r 
a d. . 

X\ = S\X\ X\ H (x 2 — Xi) 

S2 = -S2X2 + - (S m -S 2 )+- (S1-S2) 

1 — r 1 — r 

1 — a d 

x 2 = s 2 x 2 - x 2 + (xi-x 2 ) 

1 — r 1 — r 

where the output concentration s ou t is given by 

s ut = asi + (1 - a)s 2 . 

The wash-out in both compartments corresponds to the trivial equilibrium (sj„, 0, Sj n , 0), that leads to the 
trivial steady-state s* ut = Sj n . For convenience, we posit 

a 1 — a 



a>i — — , ol 2 = 



r 1 — r 

and assume, without any loss of generality that one has a 2 > a,\ (if it is not the case one can just exchange 
indexes 1 and 2). 

Remark. One has necessarily a 2 > 1 and ai < 1. 

When d — (no diffusion), the equilibrium of the system can be determined independently in the two 
compartments as simple chemostats. In this case, there is an unique globally stable equilibrium (s|,Si„ — 
s*, s 2 , Si„ — s 2 ) in the non-negative orthant, where s* = min(ai, Sj„) (i = 1, 2). 

When d > 0, we define the functions 

T 

<t>2{s\) = si + -(s in - si)(si - ax) , 
d 

1 - r 

<Ms 2 ) = s 2 H — (s m - s 2 )(s 2 - a 2 ) , 

and 

g(si) = 0i(02(si)) - Si . 

Proposition 2. When Si„ > 1 and d > 0, there exists a unique equilibrium (si,x\,s 2 ,x 2 ) of ([3]) in the 
positive orthant, where (s{, s 2 ) is the unique solution of the system 

S2 = <M s i) and 4 = 0i( s 2) > ( 4 ) 
on the domain (0, Si„) x (0, Sm), with x* = — s* (i = 1, 2). Furthermore, = = 1 when 02 = ol\ and 

ai < s* < S2 < min(a 2 , Si„) (5) 

when a 2 > oi\. 

Proof. At equilibrium, one has 

r(si + xx) + (1 - r){s 2 + x 2 ) =0, 
r(si + ii) = 0, 
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which amounts to write, from equations ([3]) 



a(si 
a(si 



x%) + (1 - a)(si 



x%) = 0, 



s\ - x\) + d(s 2 + x 2 - s\ - x*) = 0, 



or equivalently 



a 1 — a 

a + d —d 



'12 



M 



One has det(M) = a 2 — a — d < —d < and deduces the property 



S I *^']_ " — ~ "^2 — — ^ in • 

Consequently, an equilibrium in the positive orthant has to fulfill s* G [0, Sj„] for i = 1,2. Replacing by 
s j>i — s* in equations ([3]) at equilibrium, one obtains the equations 



d(s^-sl) = r(s in -s^)(s* - ax) 
d{s\~s* 2 ) = (1 -r)(s in -4)(«2 -«2) 



(6) 



which amounts to write that (s£, sj) is solution of the system ^ (see Figure 2|) or equivalently s* is a zero 
of the function g(-). 




When a 2 = ol\ = 1, one can check that s\ = s\ = \ < Sj„ is solution of (|4]). When a 2 > ai, one has 
necessarily a,\ < 1 and the condition > 1 implies g{ot\) < 0. We distinguish now two cases: 

Case a2 < Si n . If 02(^2) < Si„, notice that one has ^2(0^2) > ot 2 and then g(ot2) > 0. If 02 (0:2) > Sj n , 
notice that 4>2(ai) = a,\ < Sj n and by the Mean Value Theorem, there exists §2 G (a!i>a<2) such that 
^2(52) = Sin which implies <?(s 2 ) = s in — s 2 > 0. In both cases, one deduces by the Mean Value 
Theorem the existence of s\ £ (a%, a?) such that g(s\) = 0. 

Case a 2 > Sj ra . One has g(si n ) — with 

9 (Sin) = J5 ("1 - Sin)(Q!2 - S in ) H < . 

er d 
Rolle and Mean Value Theorems allow to conclude the existence of s± G (ai,Sj n ) such that g(s*) = 0. 
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In any case, we obtain the existence of (s*,^) solution of Q with G (at, min(a2, Sm)), that implies 
S2 = <fe(s*) > But then s\ = (f>i(s*,) < implies < min(a2,Sm)- Thus, the inequalities ([5]) arc 
fulfilled. 

Finally, notice that functions 4>i(-), 4>2{-) are both strictly concave, and steady states (s*, Sj) are intersec- 
tions of Gi, the graph of the function (f>i(-), and ^2 the symmetric of the graph of 02(0 with respect to the 
first diagonal. Consequently, if (s^s^) is a steady state different from (si„,Si„), Gi and Gi are respectively 
above and below the line segment (s*, s*,) — ( s im Sj„). We conclude that there exists at most one non-trivial 
equilibrium. □ 

Corollary 2. When Si n > 1 and d > 0, the value s± of the non trivial equilibrium is the unique zero of the 
function g(-) on (at, min(a2, s, n )). Furthermore, one has </(s*) > 0. 



Proof. When ol\ =0.2, one has s* = s| = 1 and one can easily check 

g'(st) = (1 + I( Sin - 1)) (1 + - i)) - 1 > . 

When a2 > at, one has < and we recall from the proof of former Proposition that s\ is the unique 

zero of g(-) on [at, rnin(o!2, Sjn))' We conclude that g is non decreasing at s£. Notice that 0i and 4>2 are 
concave functions and that 

MM4))) = 1 + -^-i s " L + « 2 - 2 S *) > 

implies 

We deduce that g'{si) cannot be equal to zero, and consequently one has g'(s*) > 0. □ 
The global stability of the non-trivial equilibrium is proved in the Appendix. 



Proposition 2 defines properly the map d s* ut — as* + (1 — a)s\ for the unique non-trivial steady-state, 
that we aim at studying as a function of d. Accordingly to Proposition 2, s* ut is equal to one for any value of 
the parameter d in the non-generic case a2 = a±. We shall focus on the case 012 7^ cti (and without loss of gen- 
erality we shall consider a-i > a\). We start by the two extreme situations: no diffusion and infinite diffusion. 



Lemma 1. For the non trivial equilibrium, one has 
withs?„ G (1,2). 



Proof. Under the assumptions s™ > 1 and 0:2 > at, we distinguish two cases when d = 0. 
If Sin > Qf2 ; one has s\ = ai and s\ = a-i. Then, one can write 



(1 - af (a - r) 2 



r 1 — r r(l — r) 

If Si n < a.2, on has s* = a%, s*,= Sj„ and 



> 1 



1 — eta-\ n 

'out — x ^ ^ a ws -. °in 
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(recall that assuming a 2 > «i imposes to have a < 1, and s°„ is well defined). Notice that the number s° n 
is necessarily larger than one because otx < 1, and one has also 

n (r — a) 2 

r(l-r)(l-a) 

Consequently one concludes that s* ut > 1 exactly when s^„ > s? n . Finally, remark that one has 

n r — a 2 (a — r) 2 r — a „ 

r(l — a) r(r — a) 1 — a 

Lemma 2. For s in > 1, the non trivial equilibrium fulfill 

lim sl(d) = lim s£(d) = lim s* ut (d) = 1 . 

ct— f+oo a— > + oc a—>+oo 

Proof. For any c? > 0, Proposition guarantees the existence of a unique non trivial equilibrium (s*,s£) £ 
(0, Si„) x (0, Si„) that is solution of ©. When d is arbitrary large, one obtains from © 

lim s£(d) - s* 2 (d) = . 

d— *-+oo 

From equations ((5]), one deduces also the following equality valid for any d 

r{s m - s^)(s* - ai) + (1 - r)(s in - s*)(s* - a 2 ) = , 
that can rewritten, taking into account the equality ra\ + (1 — r)a 2 = 1: 

(s m - Si)(s* - 1) = (1 - r)(s* - s*)(s m + a 2 - s* - s£) . 

Consequently, one has 

lim s*(d) = lim s^d) = 1 or lim s^(d) = lim s^d) = Sj„ . 

If a2 < Si n , the property s^ < ai valid for any c? > implies that s* cannot converges to Sj n . 

If «2 > Si n and limsj = limsj = Si„, there exists d such that rs{ + (1 — r)s^ > (sj„ + l)/2. Then, one has 

/ / * \ Kl -r) , s in + 1 - 2(rsj + (1 - r)s 2 * ) 
5 ( s i) = 75 — \ s in + on - 2s£J(s in + a 2 - 2s 2 ) H < 

GT d 

that contradicts Corollary 2. Finally, one has lims^ = lims^ = 1 and consequently lims* ut = 1. □ 

We present now our main result concerning properties of the map d i— > s* ut (d) defined at the non-trivial 
steady-state. 

Proposition 3. Assume a 2 > a%. 

- When Sm > 2, the map d H> s* ut (d) (for the non trivial equilibrium) is decreasing and s* ut (d) > 1 for 
any d > 0. 

- When Si ri < 2, the map d s* ut (d) (for the non trivial equilibrium) admits a minimum in d* < +oo, 
that is strictly less than one. Furthermore, one has 

2aia9 

d* > 



a.\ + a 2 
with s in < min(2,a2). 
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Proof. Let differentiate with respect to d the equations ^ at steady state: 
( s 2 -4) +d(d d s* 2 -d d s\) = r(s m -2s\ + a 1 )d d s\ 

" v ' 

A 

Oi - 4) + d (dd4 - ddsl) = (1 - r)(s m - 2s 2 + a 2 ) d d s* 2 



that can rewritten as follows 



A + d -d 

d -B-d 



d d s* 

ddS* 2 



(4 



Remark that one has 

A + d = d(j)' 2 {s\) 

B + d = 

det(T) = d 2 (l-<t>[(s* 2 W 2 (s* 1 )) = -d 2 g'(s* 1 ) 

From the Corollary 2, one has det(T) < and one deduces that the derivatives ddS*, ddS 2 are defined as 
follows 

I j 

d d s\ = (4-4)T7n^ 

det ( r ) (7) 



d d st 



(4-4)- 



.4 



det(T) 

Notice from inequalities (J5J) that we obtain B > and deduce dds\ > for any d. With Lemma 2 we 
conclude that Si(d) < 1 for any d. 
From equations ([?}, we can write 

o * /* M afi-(l-a)i (4 - af)r(l - r) 

OdS out = (s 2 - sj = [ai{s m - 2s 2 ) - a 2 (s m - 2s x )J - 



-dei(r) 



-det(T) 



When Si n > 2, one has A > and then ddS 2 < 0. With Lemma 2 we conclude that s£(<i) > 1 for any d. 
Then, one obtain the inequality 

o- < (s in - 2)(ori - a 2 ) < 
which proves with Lemma 2 that s* out is a decreasing function of d that converges to one. 

When Si n < 2, we write 

o- = {s in - 2)(ai - a 2 ) + 2(ai(l - s 2 ) - a 2 {\ - s*)) 

As s| and s 2 tend to one when d takes arbitrary large values, we conclude that there exists d < +oo such 
that a > for any d > d and consequently s* ut is smaller than one and increasing for d > d. We conclude 
that the map d ^ s* ut (d) admits a minimum, say at d* < +oo, that is strictly less than one. 

When d = 0, one has = a\ and s 2 = a 2 if Si„ > a 2 . Then, one obtains cr = Si n {a\ — a 2 ) < 0. So the 
map d >-> ddS* ut (d) is decreasing at d = and consequently d* > 0. 

When <i = with Si„ < a 2 , one has s 2 = Si„ and then tr = 2aia 2 — Si n {a\ + a 2 ), for which we conclude 



g < 



> 



2aia 2 

Q'l + Q'2 



Remark that this case is feasible because of the inequality 2a\a 2 < min(2, a 2 )(ai + a 2 ). We conclude that 
for Si n larger than this last value, d* is necessarily strictly positive. □ 
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Remark. The particular case a = corresponds to a configuration of a perfectly mixed tank of volume 
(1 — r)V connected to a dead-zone of volume rV. This is a way to approximate a non well-mixed tank 
or segregated bioreactors of total volume V, estimating the fraction of the volume occupied by the highly 
agitated area. 

4 Numerical computation and discussion 

Propositions 1 and 3 reveal the existence of a threshold on the value of the input concentration Si n (equal to 2 
for our choice of the parameters units) that reverses the performances of the serial and parallel configurations 
in terms of s* ut , compared to the single tank case (for which s* ut = 1): 

- for Si n > 2, there exist serial configurations such that s* ut < 1 for r large enough (i.e. the first tank 
has to be large enough), but any parallel configuration produces s* ut > 1, 

- for Si n < 2, there exists parallel configurations such that s* out < 1, while any serial configuration has 
s out > 1- There exists another threshold G (1,2) such that configurations with s* out < 1 require to 
have d large enough when Si n > sf n (cf Lemma 1). 

Furthermore, the best performance of the parallel configuration is obtained 

- for arbitrary large values of d when Si n > 2 , 

- for a finite positive d* when Sj ra G (s in ,2) (where the expression of s in is given in Proposition 3). 

For the serial configuration, the graph of the function s* ut is plotted as function of?' G [l/s m , 1] on Figure 
[5] for different values of the input concentration Si n ■ 



s 

2.0 " 


out 


Sjjj-2 


1.8 " 


S id =2.5 


\ \Sirl.75 


1.6 - 






1.4 ~ 






1.2 ~ 
1.0 






0.8 " 
0.6 





i r 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Figure 5: Comparison of s* out for the serial configuration. 

For the parallel interconnection, we depict on Figure[6]the two kind of configurations that occur, depend- 
ing on whether the number s in IS larger than one or not. 

The values of the parameters are given on the table below 





a 


r 






left figure 


0.6 


0.9 


1.14 


1.5 


right figure 


0.1 


0.9 


0.21 


1.09 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



Figure 6: Comparison of s* ut for the parallel configuration (s in > 1 on the left and s in < 1 on the right) . 

The analytic analysis of Section [3] has been conducted under the assumption of the linearity of the 
function /i(-). It is often in microbiology that the growth rate /i(-) presents a concavity, as described by the 
usual Monod (or Michaelis-Menten) function. We have computed numerically the same curves s* ut (-) than 
Figures [5] and El considering the Monod function 



65 
5 + S 



instead of the linear function (see Figure [7]). 




Figure 7: Linear and Monod growth functions. 



This function has been chosen to fulfill s* ut = 1 for the single tank configuration, guaranteeing the same 
steady state than the linear growth for this configuration. 

On Figures [8] and HI we observe that the concavity of the growth function does not change qualitatively 
the theoretical results and the existence of threshold for Si n that favourites one of the configuration. 



We notice on all the figures that the yield is better for the Monod function in the parallel configuration and 
worst for the serial one. This implies that the threshold on Si„, that was determined to be equal to 2 for the 
linear case is higher when the growth function is concave. 

Remarks. The serial configuration for the limiting value r = 1 is equivalent to a single tank. This explains 
why all the curves on Figures [5] and [8] coincide for this value of the parameter r. 
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Figure 8: Comparison of s* ut for Monod (dashed) and linear (plain) for the serial configuration. 
S ut ^ out 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1 .6 1 .8 1 2 3 4 5 6 



Figure 9: Comparison of s* ut for the parallel configuration (s in > 1 on the left and s in < 1 on the right) . 

For the parallel configuration with a = 0.1 and r = 0.9 one has ct2 = 9. This implies that for the limiting 
value d = the only equilibrium in the second tank is the wash-out when Si n < 9. This is not the case for 
the first tank but the flow rate aQ being small, the output s* ut remains closed to Si n in any case, as one can 
see on Figures [6] and [9] for small values of the parameter d. 

5 Conclusion 

Given a flow rate and the total volume of a chemostat system, this study shows the existence of a threshold 
on the value of the input concentration Si n such that above and below this threshold, serial and the parallel 
configurations are respectively the best ones with respect to the criterion of minimizing the output concen- 
tration s* ut at steady state. For the parallel scheme, the best performances are obtained for a precise value 
of the diffusion parameter that is proved to be positive when Si„ is not too small. This study concerns also 
dead-zone configurations, as particular cases of the parallel configurations. 

Whatever are the data of the problem, there always exists a configuration that is better than a single 
perfectly mixed tank. We have shown that the non-trivial steady states are unique and globally exponentially 
stable under the assumption of a linear increasing growth rate. 

Finally, this study reveals the role of the structure of the space on the performances of simple ecosystems 
or bioprocesses. The possibly non-monotonic influence of the diffusion parameter on the output steady state 
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is not intuitive, and leave further investigations open for understanding or taking benefit of this property 
for natural ecosystems (such as saturated soils or wetlands) as well as for bioprocesses (such as waste-water 
treatments) . This result can be also of interest for reverse engineering when deciding which among serial or 
parallel configurations is better ht for the modeling of chemostat-like ecosystems, providing that one has an 
estimation of the hydric capacity of the system. 
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6 Appendix: global exponential stability of the non-trivial equi- 
librium 

First, one can easily check that the domain D = fft+ is invariant by the dynamics (JXJ) and ([3]). We consider 
the 2-dimensional vector z of variables Zi = Si n — Xi — Sj 
serial and parallel configurations 



A,z 



1, 2) whose dynamics are respectively for the 



l 

1— r 



z — A.pZ — 



Notice that matrices A s and A p are Hurwitz 



-ot\ - - 
d 

l-r 



1 

l-r 



-a 2 



tr(A s ) 
ix(Ap) 



1 



r 1 — r 
-ol\ — a 2 — 



< , det(A 
d d 
r 1 — r 



d 

l-r J 



> 



r(l — r) 
< , det(Ap) = ot\a 2 



> 



So z converges exponentially toward for both systems, which implies that dynamics (TTJ) and © arc 
dissipative, in the sense that any solution of JTJ or (J3J in D converge exponentially to the compact set 
K = {{si,xi,s 2 ,x 2 ) 6 D s.t. x\ + si = Si n an d x 2 + s 2 = s m }. 

We recall a result from [19l Theorem 1.8] that shall be useful in the following. 

Theorem 4. Let $ be an asymptotically autonomous semi-flow with limit semi-flow 0, and let the orbit 
C$(t, £) have compact closure. Then the w-limit set w$(r, £) is non-empty, compact, connected, invariant 
and chain-recurrent by the semi-flow and attracts $(i,r, £) when t — > oo. 

6.1 The serial configuration 

Proposition 5. Under the condition s; Tl > 1/r, any trajectory of (TTJ with initial condition in D such that 
(si(0),xi(0)) 7^ (si„,0) converges exponentially to the unique non-trivial steady-state (s*, x*, s 2 , x 2 ) given 
by Proposition 1. 

Proof. Dynamics (fTJ) has a cascade structure. It is straightforward to check that the solutions of the 
(si,xi) sub-system converges asymptotically towards the non-trivial equilibrium (1/r, s; n — 1/r) from any 
initial condition away from the wash-out equilibrium (si„,0). From the convergence of z 2 toward 0, we 
deduce that the s 2 variable has to converge to the bounded interval [0, Si n ] and that its dynamics can be 
written as a scalar non autonomous differential equation: 



S2 



"S 2 (Sm -8 2 - Z 2 (t)) + 



(8) 
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This last dynamics has the property to be asymptotically autonomous with the limiting differential equation: 

1 



h = f(s2) = -s 2 {s in - s 2 ) + 



1 - r 



(1 A - «a) 



(9) 



Statement of Proposition 1 implies that this last scalar dynamics has a unique equilibrium s 2 that belongs to 
[0, Sin]- Furthermore, one has /(0) > and f(si n ) < 0. Consequently any solution of (|9|) in [0, Si n ] converges 
asymptotically to s 2 . Then applying Theorem 4, we conclude that any bounded solution of dHJ converges to 
s 2 . Finally any solutions of the (s 2 ,x 2 ) sub-system converges asymptotically to (s 2 , s,„ — s 2 ). 

The Jacobian matrix of dynamics |T]) at the non-trivial equilibrium (s*, x\, s 2l x 2 ) is of the following form 
in {zi, z 2 , s\, s 2 ) coordinates 



A s 
★ J* 



with J* 



l 

l-r 







2s% — rp— 



Recall that A s is Hurwitz. The eigenvalues of J* are — (sj n — s*), — (sj n — s 2 ) — (1/(1 — r) — s 2 ) that are both 
negative numbers, accordingly to Proposition 1. The exponential stability of the non-trivial equilibrium is 
thus proved. □ 

6.2 The parallel configuration 

Proposition 6. When s;„ > 1 and d > 0, any trajectory of (J3)) with initial condition in D such that 
xi(0) > and x 2 (0) > converges exponentially to the unique non-trivial steady-state (s*, x\, s 2 , x 2 ) given 
by Proposition 2. 

Proof. Considering the time vector z(-), the (si, s 2 ) sub-system of dynamics © can be written as solution 
of a non-autonomous planar dynamics 



si 
si 



Sl{zi(t) + Si - s in ) + ati(s in - Si) + ^(s 2 - Si) 
s 2 (z 2 (t) + s 2 - s in ) + a 2 (s m - s 2 ) + jz^(si - s 2 ) 



(10) 



We know that z converges to and consequently the vector S of variables sj., s 2 converges to the set 
S = [0, Si n } x [0, Si n ]. We study now the limiting autonomous dynamics 



si 



si)(a 1 - si) + ^(s 2 ~ s x ) 
s 2 )(a 2 - s 2 ) + j4^(si - s 2 ) 



(11) 



on the domain S. Let B be the boundary {s\ — Si n } U {s 2 = Sj n }. On the domain S \ 23, we consider the 
vector a of variables <n — log(st„ — Si), whose dynamics can be written as follows 



a = F(a) = 



—ai + Si 

-a 2 + s in 



e ' 

^2 



(1 



1— r 



(12) 



One can easily compute 



div(F) 



-e" 
r 



1 



< 



From Poincare-Bendixon theorem and Dulac criterion, we conclude that bounded trajectories of (|12|) cannot 
have limit cycle or closed path and necessarily converge to an equilibrium point. Consequently, any trajectory 
of (fTTj) in S either converges to the rest point S* = (s*, s 2 ) or approaches the boundary 23. Notice that one 
has 



^ini &j ^ $in ~^ $i ^ 



(i ¥> J) 
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So the only possibility for approaching B is to converge to the other rest point S° — (sj n , Sj n ). This shows 
that the only non-empty, closed, connected, invariant and chain recurrent subsets of S are the singletons 
{S*} and {S }. 

Applying Theorem 4 we conclude that any trajectory of (jlOp . issued from initial condition of dynamics 
in T>, converges asymptotically to S* or S°. Consider now any initial condition with xi(0) > and 

£2(0) > 0. We show that the solution (si(-), S2O) of ([3]) cannot converge to 5°. If it is the case, there exists 

T < +00 such that one has 

S\(t) > a\ and rsi(t) + (1 - r)s 2 {t) > 1 for any t > T 
under the assumption Si n > 1. Let us consider the function 

V(xi,x 2 ) = min(rxi + (1 - r)x%,xx) 
(see Figure [T0|) and v(t) = V(xx(t),X2(t)) that is positive and tends to when t tends to +00 



Figure 10: Iso- value of the V function. 

If X\{t) < X2(t), one has v(t) = x\(t) and 

v = x\ > (si(t) - a\)xi > for t > T 

If x\(t) > X2(t), one has v(t) — rx\{t) + (1 — r)x2(t) and 

v = rxi + (1 - r)x 2 = r(si - ot\)xi + (1 - r)(s 2 - a 2 )x 2 
> (rs 1 + (1 - r)s 2 - l)x 2 > for t > T 

We conclude that the function t 1— > is non-decreasing for t > T and consequently cannot converge to 
zero, thus a contradiction. 



The Jacobian matrix of dynamics ([3]) at the non-trivial equilibrium (s*, S2, ^2) i s 01 the following form 
in [z\, Z2, si, S2) coordinates 



A P 
★ J* 



Recall that A p is Hurwitz. One has 



with J* — 



dfit(^) = ^j^j(0i(*a)^(*i)-l) and tr(J*) = -^ 2 (^) - ^^(4) 
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The function ^i(-) being concave, one has </>i(si n ) < 0i(s£) + (^(sJXsm — sj)- Along with the inequalities 
s in > s\ and (f>i(si n ) — s.;„ > s* = ^(s^), one deduces (^[(st,) > 0. Recall from Corollary 2 that one has 
9'( s i) = fiii^^i 8 !) — 1 > 0. Then the inequality ^(s*) > is necessarily satisfied. Finally, we have 
shown det(J*) > and tr(J*) < 0, that guarantee the exponential stability of the non-trivial equilibrium 

Remark. The wash-out equilibrium (sj n , 0, Sj n , 0) is not necessarily hyperbolic. This explains why we cannot 
use the Convergence Theorem for asymptotically autonomous dynamics given in Appendix F of |29j . 
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